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1 .0 INTRODUCTION 


This report is a summary of the past year’s effort under 
Contract NAS8-33698 , entitled "Transient and Diffusion Analy- 
sis of HgCdTe" . The primary purpose of this work is to provide 
analysis and computational support for a proposed Materials 
Processing in Space flight project involving crystal growth 
of mercury cadmium tell u ride alloys under microgravity condi- 
tions. The overall goal of this flight project is** to obtain 
homogeneous material with minimum radial and axial composition- 
al variations. Accordingly, this work addresses the problem 
of solute redistribution during directional solidification of 
HgCdTe . 

This work reports on new results obtained using the one- 
dimensional model, the interface deconvolution model, and add- 
resses the problem of two dimensional diffusion and shape in- 
stability. 

Two papers have been submitted for publication and pre- 
liminary copies of these will be included in the appendix. 

In addition, computer code for the models discussed here is 
documented and issued under a separate cover to the appropriate 
NASA personnel. This report contains experimental data from the 
work of Dr. D.C. Gillies, formerly with USRA/MSFC , and is in- 
cluded for comparison with the model predictions. 


Page 1 



2.0 BACKGROUND 

The solid solution semiconductor HgCdTe has emerged as 
the premier variable energy gap intrinsic infrared detector 
material for the future. Considering Hg 1 Cd Te as a mixed 
crystal of HgTe and CdTe, the material properties, including 
the bandgap, may be varied by changing the solute (CdTe) mole 
fraction x. The bandgap varies linearly with x from approxi- 
mately -.3eV (x=0 , semimetal) to 1 . 6eV (x=l, CdTe) allowing 
selection of a suitable proportions of HgTe and CdTe. 

f 

It is this property, however, that is the source of one 
of the primary problems in the growth of high quality HgCdTe. 
Requirements for large area focal plane arrays for IR detectors 
necessitate compositional uniformities of Ax< . 0 1 or less over 
a crystal slice of diameter lcm. or larger. Compositional 
variations, either radial or longitudinal, lead to variations 
in bandgap and thus detector response over the area of the 
detector. Growth techniques that eliminate or minimize com- 
positional variations are of primary importance in producing 
state-of-the-art material capable of meeting detector performance 
requirements . 

In order to produce more homogeneous material we must first 
understand and quantify the causes of compositional variations in 
HgCdTe. Although several growth techniques have been invest- 
igated by other workers, we concentrate on growth by directional 
solidification as this offers a high potential yield of material 
from a relatively simple growth technique. Directional solidifi- 
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cation of HgCdTe results in a solute (CdTe) profile that is 
characterized by a mean axial solute distribution with an 
initial transient region, a steady-state or near steady- 
state region depending on ampoule pull rate, and a final 
transient region characterized by a decreasing CdTe con- 
centration as the end of the ampoule is reached. This 
indicates that the primary means of solute redistribution 
in HgCdTe is via diffusion through a diffusion layer ahead 
of the solid-liquid interface. 

t 

Concentrating on the axial solute profile, SEMTEC has 
developed a one-dimensional model based on the diffusion 
equation to predict the axial solute distribution in a 
finite length ampoule. The model assumes a planar inter- 
face and linear thermal field, but accounts for variations 
in the segregation coefficient and interface temperature 
with composition, factors that are particularly important 
in HgCdTe. The model also calculates a growth velocity 
determined by the thermal field and interface concentration 
rather than assuming a constant growth velocity equal to 
the pull rate. 

The inputs to the one dimensional model are the phase 
diagram (pseudo-binary) , ampoule pull rate and length, initial 
solute concentration, and the effective diffusion constant 
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(taken independently' of temperature and concentration). The 

model gives axial solute concentration in the solid, interface 

position, and interface velocity at each chosen time interval. 

By comparing experimental solute profiles to the model results 

for various values of the diffusion constant, we have arrived 

-52 -52 

at a value of D between 5x10 cm /sec and 7x10 cm /sec for a 
wide range of ampoule pull rates (.076 cm/hr to .41 cm/hr). 

The degree of agreement between experiment and theory is quite 

good, confirming that the axial solute distribution is dominated 

/ 

by diffusion for HgCdTe. 

Details of the model are given in a paper submitted for 
publication and included in the appendix of this report. Results 
obtained this year using the model will be addressed in a later 
section of this report. 

Further examination of the solute distribution in Bridgman 
growth HgCdTe reveals that there is a significant radial variation 
in the solute. This radial variation is present in crystals 
grown under a variety of thermal conditions and growth rates. 
Although axial solute variations are an unavoidable consequence 
of diffusion limited growth, it should be possible to eliminate 
radial variations under the proper. growth conditions, at least 
during steady state growth. The causes of this radial solute 
variation is a key issue in the ground based HgCdTe research 
effort . 

Experimental results to date indicate that in the steady 
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state growth region, the isoeoncent rat ion lines are concave 
(concentration of CdTe at r=o is less than at the crystal outer 
surface) . The isoconcentration lines are not perfectly axisymmet- 
ric, indicating mis alignment of the ampoule axis with the g-vector 
and/or thermal field. 

It is appropriate at this point to indicate some fundamental 
observations regarding solute redistribution. In the steady state 
growth region we have (by dif inition) no variation in the axial 
solid composition at a fixed radius, i.e. the solid' concentration 
is a function of the radial variable r alone. If the interface is 
an isotherm, then it is also an isoconcentrate and homogeneous 
material results. Conversely, if the interface is found to be an 
isoconcentrate, then it must be an isotherm and homogeneous mater- 
ial results. As a corollary, if radial composition variations are 
present in the steady-state region, the interface must be neither 
an isotherm nor an iso concentrate, a situation which complicates 
the analysis of the growth behavior. 

Since radial compositional variations are observed in HgCdTe 
samples, we must seek the cause of these variations in order 
to ascertain whether appropriate growth techniques may eliminate 
them. Several factors may influence the steady-state solute 
distribution; the thermal field, interface shape, convect- 
ion, and diffusion may all play a part in determining steady 
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state conditions. An additional mechanism for solute redis- 
tribution in the HgCdTe system has been postulated. Consider 
a system such as HgCdTe in which the rejected component (HgTe) 
is more dense than the bulk fluid. As a diffusion layer is 
built up during growth, a layer of HgTe rich liquid is estab- 
lished just ahead of the solid-liquid interface. If the inter- 
face becomes curved, say concave, the heavier HgTe rich liquid 
will flow along the interface and concentrate at the lowest 

point along the interface, leaving the center of the interface 

( 

HgTe rich as compared to the outer edges. As HgTe is also the 
lower melting point element in the pseudobinary system HgTe- 
CdTe , the solidus temperature of the center of the interface is 
lowered with respect to the interface edges (for a concave 
interface). As the temperature increases as we move away from 
the interface into the liquid, this solute profile will result 
in an increase in curvature (more concave) in the interface. 

This is a dynamic process that once begun leads to interface of 
increasing curvature limited only by thermal consideration , 
Similar arguments apply if we consider a convex interface with 
accumulation of HgTe at the edges of the sample. This "shape 
instability" was first postulated for the HgCdTe system by David- 
son (1) as a possible cause of radial compositional . variations , 

We will briefly address each of the above mentioned mech- 
anisms for radial solute redistribution in the following pages 
and provide model results and experimental details in later 
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sections of this report. 

The thermal field is important initially as it determines 
the location and shape of the initial interface. Assuming 
nucleation and faceting effects do not effect the growth pro- 
cess after some time interval, then the first to freeze portion 
freezes out of a liquid of composition near C q (taken to be x = 
.20 in most experiments) at a temperature of nearly 800°C. 
Furnace temperatures are generally chosen to yield the steady 
state solidification isotherms (705°C for C o =.2Q) at the middle 
of the adiabatic zone, and thus one might initially suspect 
that the first to freeze portion occurs above the adiabatic zone 
with a convex isotherm. However, initially the thermal coupling 
between the ampoule and the cold zone in the furnace is poor 
since only a small portion, if any, of the ampoule is in the 
cold zone. This shifts the isotherms downward in the furnace 
and may result in initial freezing well down in the cold zone. 
Thermal problems of this nature are difficult to model since 
the heat transfer characteristics are highly dependent on how 
the ampoule is constructed and supported in the furnace; ex- 
periments using model materials in an instrumented ampoule are 
needed to determine the early thermal behavior. 

Independent of the location in the furnace that the first 
to freeze occurs, it is evident that it is difficult to set 
experimental conditions such that the solid-liquid interface 
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remains on a planar isotherm. Given that the interface likely 
becomes curved on approach to steady-state , shape instability 
and/or lateral diffusion along a curved interface will result 
in radial compositional variations in the steady state region. 
Due to the high vapor pressures encountered in HgCdTe, it is 
difficult to instrument the ampoule with thermocouple or elect- 
rodes for Peltier pulsing. We must use other techniques to 
determine the interface shape and concentration. 

We find ourselves in the position of realizing that the 
steady-state solute distribution depends greatly on the initial 
stages of growth, and yet cannot adequately predict the initial 
behavior due to complicated thermal fields and nucleat ion be- 
havior. Thermal modeling of the steady-state region in which 
isothermal boundary condi tons can be applied to the ampoule 
some distance away from the interface is on a firmer footing. 

We also recognise that chemical analysis indicates that there 
are radial composition variations in the steady-state, thus the 
interface is neither given by the shape of the isotherms nor the 
isoconcentrates. Since we may predict with some confidence the 
thermal field at steady state and experimentally measure the 
concentration field in the grown crystal, we may use this in- 
formation to "deconvolu te ,r the existing interface shape at any 
point in the steady state region. 
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Details of the deconvolution technique are given in 
a paper submitted for publication which is included as an 
appendix to this report. Examples of the results obtained 
from the deconvolution process are given in a later section 
of this report. The basic principle of interface decon- 
volution is based on the observation that the solid-liquid 
interface is the locus of all points in the sample whose 
solidification temperature coincides with the thermal field 

in the sample at a particular time. Let C (x,y) be the solid 

s 

composition profile measured on the crystal. The solidus 

curve of the phase diagram may be used to obtain (S(C (x , y ) ) 

s 

where S is the solidus temperature corresponding to the con- 
centration C (x,y) and we assume no undercooling. By approp- 
s 

riate use of a thermal model we may calculate the thermal 

field present in the crystal T(x,y), We thus have two tempera 

ture distributions in the crystal at any point in time, the 

temperature at which each point froze S (G (x , y ) ) and the cal- 

s 

culated temperature at each point in the crystal T(x,y). The 

interface is determined as those points at which S (C (x , y ) ) is 

s 

equal to T(x,y) . 

In actual practice we must use an iterative procedure 
to accurately determine the interface location since the inter 
face represents a change in thermal conductivity and is thus a 
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boundary condition to be applied in the thermal model. Details 
of this iterative procedure are given in the appendix. The 
result of the deconvolution process is an interface location 
and solute concentration along the interface that: is : self- 
consistent with the thermal boundary conditions and measured 
solute concentration profiles. 

Once a steady-state interface shape and composition have 
been established, we may seek to discover the causes of the 
compositional variations along the interface. It is known 
that lateral diffusion along a curved interface may lead to 
radial variations in solute redistribution. Coriell, Sekerka, 
et al (2) have addressed the problem of solute redistribution 
by diffusion along a curved interface. They assumed an inter- 
face shape and calculated the resulting interface solute pro- 
file due to lateral diffusion. Note that it is not necessary 
to determine the thermal field to solve this problem; the 
assumed interface shape and resulting composition may not be 
consistent with thermal boundary conditions encountered in 
an experiment. The work of Coriell and Sekerka has been applied 
to HgCdTe and sample calculations are given in a later section 
of this report. Based on preliminary experimental work and apply- 
ing Coriell and Sekerka f s calculations, it does not appear that 
lateral diffusion alone can account for the magnitude of radial 
segregation found in HgCdTe samples. -We must then address the 
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mechanism of shape instability to ascertain if gravity driven 
flow along the interface is the cause of these compositional 
variations. In a later section, the appropriate equations are 
developed and certain conservation conditions are derived 
pertaining to shape instability. 

The work thus far has followed a logical sequence. The 
one dimensional model is applied to determine an effective 
diffusion constant and to investigate the influence of various 
processing parameters such as pull rate and gradient on the 
axial profile. Two dimensional interface deconvolution is used 
to estimate the magnitude of interface curvature that leads 
to radial inhomogeneities. Once an approximate interface shape 
is established, we can apply two dimensional diffusion calcul- 
ations to determine if diffusion is the dominant mechanism for 
radial solute variations. Finally we address the problem of 
shape instability and develop appropriate equations suitable 
for numerical solution. In the following sections we present 
calculat ional and experimental details concerning each of these 
topics. 
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3.0 ONE DIMENSIONAL TRANSIENT ANALYSIS 

Details of the one dimensional transient analysis program 
are given in the appendix to this report. Previous results 
are also given in a report for contract NAS8-33698 dated 
20 December 1980. 

3 . 1 RESULTS USING ONE DIMENSIONAL TRANSIENT ANALYSIS 

The first objective of the 1-D program was to determine 
the effective diffusion coefficient. Per this the original 
intention was to utilize both initial and final transient re- 
gions of crystals grown at constant translation rate. For 
reasons to be explained this procedure did not lead to good 
results with the initial transient and other procedures in- 
volving changes in the translation rate and stopping the crystal 
during growth were devised. Some of the most recent uses of 
the program will be described in the following sections. The 
implication of the results on the design of future growth ex- 
periments will be stressed. 

3.2. UNINTERRUPTED DIRECTIONAL SOLIDIFICATION 

Many sets of data have been calculated using as variables 

the translation rate, the length of the ampoule, the effective 

diffusion coefficient and the nominal composition. 

Results typical of what the model can achieve are shown 

in Figures 1-4. To illustrate the program a 20cm crystal of 

Hg n Cd Te was "solidified" in a gradient of 200 °C /cm at 
u.oUU.20 

a constant translation rate of 2mm/hour . Figure 1 shows the 
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Figure 1 
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composition along the length of the crystal assuming diffusion 

-52 -52 

coefficients of 5x10 cm /sec and 7x10 cm /sec, A clear dis- 
tinction can be seen particularly in the region of the final 
transient. The value during the so-called steady state region 
never quite reaches the nominal value of 0.20; this is due to 
series truncation which builds up during the numerical analysis. 
Absolute accuracy must be traded off against excessive program 
running time to reach a r e p r e s e n t a t ive result . 

Other data from this model are shown in Figure 2. The 

i 

location of the interface demonstrates solidification initiating 
high up in the hot zone of the furnace and proceeding down to the 
mid-plane (0 on the graph) where steady state sol id j f ica t i on takes 
place. In the final transient the interface moves down into the 
cold end where the mercury rich material freezes. 

Figure 2 also illustrates the interface velocity which reach- 
es a maximum equal to the translation rate only during steady state 
solidification and during the very last moment of crystal growth 
when the melt is completely depleted of cadmium. The decrease 
and subsequent increase in the velocity during the final trans- 
ient occurs as a result of the inflexion point in the composition 
data as shown in Figure 1. As the final transient regime is 
entered the interface moves down into the cold zone to solidify 
liquid which can no longer be replenished by diffusion with cadmium. 
The interface thus moves slower than the pull rate. However towards 
the end of growth the liquid composition tends asymptotically 
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to the mercury tellur ide composition which will grow at the 
translation rate. The interface velocity thus increases. 

The effect of varying the temperature gradient was also 

calculated. Several different, but plausible, temperature 

gradients were input as data and the resulting transient analysis 

curves for 100 °C /cm and 400 ° C / cm are shown in Figure 3 for a 

sample pull rate of 2mm/hr . and a diffusion constant of 7x10 
2 

cm /sec. In comparison to varying the diffusion coefficient, 
which has a large effect on the results, it can be seen that the 

i ■ 

shapes of the transient analysis curves are not altered much by 
gradient. This result has important implications, in that the 
gradient within a silica ampoule loaded with HgCdTe cannot be 
measured with any degree of accuracy. Thus, a reasonable value 
for the effective diffusion coefficient can be obtained due to 
the high sensitivity of D and the low sensitivity of G, the 
temperature gradient. 

There is one other variable which enters into the transient 
analysis equations, namely translation rate. Figure 4 shows the 
effect of varying the translation rate over several orders of 
magnitude. These curves are important in that they give immediate 
information as to what the minimum translation rate must be to 
expect any steady state growth in the crystal . There are also 
implications here as to the desirable length of ampoule and the 
trade off between translation rate and ampoule length. A trade 
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off exists between translation rate, thermal gradient, and ampoule 
length. As the translation rate is increased, the thermal gradient 
necessary to avoid constitutional supercooling increases. On the 
other hand, decreasing the pull rate results in longer transient 
regions with a c o r r e sp ond ing shorter steady state region for a 
fixed ampoule length. The one dimensional model is well suited 
for choosing proper growth parameters prior to actual growth 
experiments . 

Data pertaining to a real crystal together with the compos- 
itional data of the same crystal are shown in Figure 5. This 
crystal was grown at NASA by Dr. L.R. Holland of Athens State 
College. Translation rate of 0.762mm/hr is slower than the 2 nun/ 

hr used in the example of Figure 1 and the transient regions are 

-5 2 ' 

somewhat longer. For the diffusion coefficient of 7x10 cm /sec 
it is evident that there is little steady state region* The 
crystal exhibits end effects before the initial transient is 
completed. Assuming this diffusion coefficient, a much longer 
crlystal than this would have to be grown to ensure a reasonable 
length grown in the steady state regime. Compositional data for 
the real crystal demonstrate a reasonable fit except in the init- 
ial transient region. There are several reasons for the initial 
fit being bad. The two main ones are as follows. 


1 . The 

ampoule 

was designed 

with a 

conical end in 

order 

to 

encourage the 

growth 

of a single 

crystal 

from the first 

to f r 

eeze 

nucleus . The 

mode], , 

however , ass 

ume s a 

flat interface 

and in 

the 
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initial transient region during which the interface moves towards 
the cold end there is little chance of this occurring especially 
with a changing cross-section* 

i 

2. It was evident from other results that nucleation did 
not take place at the tip of the ampoule. The starting point for 
unidirectional solidification is thus ill defined, making compar- 
ison to model results difficult* 

Once nucleation has taken place and solidification occurs 

in the cylindrical section of the ampoule the experimental data 

fall much closer to the model. The experimental results shown 

here represent a mean of nine points taken with a scanning electron 

microscope fitted with energy dispersive spectrometry. As these 

readings are taken on a cross-section a flat interface is assumed. 

However, radial solute variations do exist and thus these data 

represent the mean axial solute distribution. The final transient 

-5 2 

shows a good fit to the diffusion coefficient of 7x10 cm /sec . 
Further details of this fitting may be found in the paper included 
in the appendix. 

The use of the model in the direct mode results in determining 
effective diffusion coefficients, optimum translation rates and 
crystal lengths for obtaining a reasonable length of crystal grown 
in the steady state. An assumption is made here as to a known ther 
mal gradient and, more important, a flat interface. For a given 
gradient it is possible to approximate the fastest translation 
rate before breakdown takes place. By applying the classical con- 
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stitutional supercooling criteria, we find a temperature gradient 

of 1 7 0 ° C / c m is needed to grow cry st als at 5mm/hour assuming D=7xl0 
-5 2 

cm /sec. It was found, however, that a crystal grown at 4mm/hour 
exhibited dendritic growth. There may have been several reasons 
for this. A gradient of 1 8 0 ° C / cm was measured in an empty silica 
ampoule. In reality the gradient could be a lot less; in situ 

instrumentation of mercury cadmium tellur ide ampoules has never 
been achieved. Breakdown was observed in the initial transient 
region and it is unlikely that the interface was flat. More 
details of this crystal are given in Section 3.3* 

3 . 3 DI RECTIONAL SOLIDIFICATION WITH R A TE CHANGES 

The transient analysis model can be applied to simulate a 

change in translation rate during growth. Assuming a temperature 

- 5 2 

gradient of 1 8 0 ° C / c m and D~7xl0 cm /sec several rate changes 
were modeled. The translation rate of 4mm/hour was thought to 
be the fastest that could be used without breakdown and this was 
used in one simulation with a rate change to 1mm /hour after steady 
state had been reached. The data obtained are shown by the solid 
line in Figure 6, Following the rate change a new diffusion layer 
is set up consistent with the new translation rate. In the ex- 
ample shown the crystal is too short for growth to even approach 
steady state, end effects are observed immediately, and the com- 
position varies continuously. 

The simulation above models a real, crystal grown under these 
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conditions by Dr. F.R. Szofran at MDRL. The mean analysis at 
various locations along the length of the crystal' is also shown 
in Figure 6. As explained in Section 3.1 a translation rate of 
4mm/hour led to breakdown of the interface with a resulting 
dendritic structure. Concentration values varied considerably 
across a diameter and there is no guarantee that a valid mean 
value can be determined. 

Nevertheless the analysis points from the real crystal do 
bear some resemblance to the model. Following the rate change, 
growth of single crystal material took place. The existence of 
a second initial transient is clear. As this material would be 
prenucleated such a region should give accurate diffusion data 
provided the initial growth had a more or less flat interface. 

Such an experiment is thus worth repeating providing the initial 
growth takes place under more controlled conditions. 

This experiment has other implications. The determination 
of the shape and location of the interface is of paramount im- 
portance as a first stag e of designing a controlled growth ex- 
periment. Due to the impracticability of instrumenting a mercury 
cadmium telluride ampoule indirect methods have to be used. The 
deconvolution method discussed earlier is one such method' another 
is the use of interface perturbation such as is present in a rate 
change. In principle it should be possible to determine the 
location in the crystal of the interface at the time the rate 
change occurred. At such a position there is a discontinuity in 
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the chemical composition caused by the setting up of a new 
diffusion layer. In the case of crystal CMT-1 shown in Figure 
6 this determination was not attempted due to the breakdown of 
the interface caused by constitutional supercooling. This 
technique should be further investigated under more favorable 
experimental conditions. 

3. 4 INTERRUPTED GROWTH 

The model can also be used to simulate behavior when trans- 

i 

lation of a crystal is stopped entirely. Under these circumstances 
two main changes occur assuming diffusion in the solid state can 
be neglected. First the liquid will homogenize reaching a compos- 
ition somewhat lower in cadmium telluride than the original compos- 
ition. As this liquid will no longer be in equilibrium with the 
solid, growth will continue and the interface will advance (to- 
ward the hot end). During this period the solid being deposited 
will become progressively richer in cadmium telluride until it 
is in equilibrium with the liquid, now homogeneous. At this moment 
the interface will lie along an isotherm within the furnace and 
an ideal situation for growth from a prenucleated crystal exists. 

The model is useful here in demonstrating the lengths of ampoule 
needed to perform such an experiment and still achieve steady state 
following a re-start of growth. The length of the interruption can 
also be determined. 
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An example illustrating the method and the effect of two 
different stoppage times is shown in Figure 7. In this case a 
rate change was incorporated as well. It should be noted that 
the steady state region following the restart lies at a lower 
value than the first steady stated this is due to the incorpor- 
ation of cadmium rich material in the initial transient. It is 
strongly recommended that a controlled experiment such as this 
be attempted as soon as possible. Once again, interruption 
of growth with or without a rate change offers a means of locating 
the interface due to the change in composition whern growth is 
re-initiated. 
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4.0 INTERFACE DECONVOLUTION 


4 * 1 BACKGROUND 

As described in Section 2, interface evaluation is a 
critical factor in Bridgman St ockbar ger crystal growth. The 
location, shape, composition, temperature and movement of the 
interface during the experiment are important parameters, crucial 
for the control and design of a successful crystal growth program. 

In this context mercury cadmium tell u ride poses its own problems. 

The mercury vapor pressure at the temperatures under consideration 

f " 

is in excess of 20 atmospheres and it has so far proved impossible 
to contain the material in an ampoule instrumented with feedthrough 
leads . 

Thus the technique of Peltier pulsing which has so success- 
fully been applied to doped germanium, and is at present being 
developed for lead tin telluride, cannot be used here. Any infor- 
mation about the interface has to be gathered either by more 
complex analytical techniques, such as the deconvolution method 
to be described, or by experimental perturbations such as the 
rate change described in Section 3.3. Any perturbation technique 

will lead only to one experimental data set for interface character- 
ization, Any movement or change of shape of the interface during the 
crystal growth is difficult to ascertain in a single experiment. 

This section reviews the work done so far in developing the 
deconvolution model, describes typical results and presents ideas 
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suitable for modifying the model to assist in the planning of 
future experiments* In another sub-section results achieved 
with the perturbation techniques are described and an assess- 
ment is made of their utility in future experiments. 

The deconvolution technique was originally conceived by 
Dr. D.C* Gillies at NASA, while working as a USRA Visiting 
Scientist. SEMTEC has formalized the concept into a mathematically 
rigorous form and developed a computer program that performs the 
deconvolution process. The iteration technique described in the 
paper in the appendix is necessary to give accurate results, but 
the general characteristics of the process may be described below. 

4 . 2 DECONVOLUTION - APPLICATION 

The deconvolution technique effectively determines the prior 
locations of the liquid-solid boundary after growth is completed. 
For a mathematical description of the methods the reader should 
refer to the appendix in this report-. The first attempts at separa- 
ting solidus temperatures from thermal field were made by hand 
using transparent paper. They serve to illustrate the method and 
how the computer techniques evolved from it. 

A schematic representation of the procedure is shown in 
Figure 8a-c. Figure b is the assumed thermal field which in this 
case is decoupled from the crystal. The midplane is at 705°C, the 

solidus temperature for Hg n on Cd on Te * In Figure 8a are shown 

U • o U U * U 
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isoconcentration lines in the grown crystal together with their 
equilibrium solidus temperatures. These curves have to be con- 
structed from the analysis data of the crystal. The super imposition 
of these two graphs demonstrates within the furnace frame and the 
crystal frame locations of the solidification front. For one in- 
stant this procedure is shown in Figure 8c. The coincidence of 
thermal field and solidus temperature is shown at the circles. 

With additional isotherms the complete interface can be reconstructed 
as demonstrated by the crosses. It can readily be seen that the 

r 

interface is not along an isotherm nor is it an iso concentrate. 
Furthermore, in this case, the interface reflects the curvature 
of the solidus and not that of the thermal field. The movement of 
Figure 8a and 8b relative to each other is a function of time and 
the translation rate of the ampoule or crystal. At any time a 
different set of coincidence points can be found corresponding to 
that instantaneous location of the interface,.. The above procedure 
lends itself to a computing technique; the program and its use 
are included as an appendix. 

4 - 3 DECONVO LUT I ON - EXAMPLE 

As explained previously, the complete interface deconvolution 
process requires iteration between the thermal model and decon- 
volution procedure to obtain an interface location and composition 
that is consistent with both the thermal boundary conditions and 
the measured solute profile. We may illustrate the deconvolution 
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procedure using composition data from an actual crystal and apply 
various assumed thermal fields to qualitatively describe the inter- 
face behavior . 

The crystal is that used to illustrate the 1-D transient 

analysis model in Section 3.2. Crystal composition was Hg n on 

U . o U 

C d o 2 Q Te » t * ie diameter was 8mm, the length 21.2cm and the trans- 
lation rate was 0 . 7 6 2mm/ ho ur . The complete crystal analysis 

C ( x , y ) was obtained by energy dispersive spectrometry. The raw 

■/ 

data for selected radial values of x were fitted to a polynomial, 
as were the thermal data, in order to implement the deconvolution 
procedure. C(x,y) values were converted to S(x,y), the correspond^ 
ing equilibrium solidus temperatures by the deconvolution computer 
program. It is S(x,y) that has been plotted in Figure 9 to ill- 
ustrate compositional variations throughout the crystal. This 
crystal was not grown under ideal conditions; there is a pronounced 
asymmetry and a pronounced drop in the isoconcentration lines close 
to the right hand edge. 

Of significance is that the iso concentration lines are ex- 
tremely curved even towards the top of the illustrated section 
which lies close to the steady state region. If the interface 
were an isoconcentrate its curvature would have to be extreme. 

It is conceivable that the thermal field is distorted to the ex- 
tent that solidification along an isotherm could yield the iso- 
concentrate shapes presented in Figure 9; however, we will see 
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that postulating a more reasonable thermal field leads to an 
interface of moderate curvature consistent with the composition 
data. 

Two thermal fields were chosen to illustrate the application 
of the method. The two fields are shown in Figure 10. Each has 
one flat isotherm and a temperature gradient of 200°C/cm along 
the centerline of the furnace. In one case the flat isotherm 
is at 700°C and in the other at 750°C. These two values were 
chosen as extremes. In thermal field 1 (700°C flat isotherm) 
solidification would take place above the flat isotherm in a 
region where the isotherms are convex (i.e. directed upward at 
the center). Conversely in thermal field 2 the isotherms are in 
the opposite direction. 

The "deconvolution 11 will give resulting interface shapes, 
locations and temperature anywhere within the range of the cal- 
culated polynomials for S(C(x,y)) and T (x , y ) . Figure 11 shows 
the interface for both thermal fields at two selected points in 
the crystal and thus at two different times during growth, measured 
in the furnace frame. The temperatures and actual locations are 
given in Table 1. The interface at the same time for thermal field 
2 is also shown in Figure 11 and in Table 2. It can be seen that 
in both cases the interface was concave (i.e. directed downward at 
the center), but the degree of curvature is very dependent on the 
isotherm curvature within the region of solidification. Thermal 
field 1 has convex isotherms at the solidification temperature 
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TABLE 1 


* 

Time 
1 Hour 


Time 
10 Hours 


Time 

100 Hours 


RESULTS OF DECONVOLUTION USING THERMAL FIELD 1 


x ,mm 

Temperature , °C 

y ,mm 

y ,mm 

-3.916 

721.9 

26.01 

0.96 

-3.440 

720.6 

25.97 

0.93 

-2.919 

718.6 

25.90 

0.86 

-2.252 

716.1 

25.82 

0.78 

0 

714.2 

25.75 

0.71 

2.252 

714,6 

25.75 

0.71 

2.919 

714.8 

25.72 

0.69 

3.440 

715,8 

25.76 

0.72 

3.916 

718.5 

25.86 

0.82 

-3.916 

719.1 

32.72 

0.84 

-3.440 

717.7 

32.68 

0.80 

-2.919 

716.1 

32.63 

6.74 

-2.252 

713.9 

32.55 

0.67 

0 

712.4 

32.02 

0.62 

2.252 

712.2 

32.47 

0.59 

2.919 

713.3 

32.50 

0.61 

3.440 

713.6 

32.50 

0.62 

3.916 

715,8 

32,58 

0.70 

-3.916 

711.6 

100 . 79 

0.51 

-3.440 

710.2 

100.75 

0.46 

-2,919 

708.9 

100.69 

0.41 

-2.252 

707.8 

100.66 

0.38 

0 

706.4 

100.60 

0.32 

2.252 

707.8 

100.66- 

0.38 

2.919 

706.6 

100.58 

0.30 

3.440 

708.6 

100.67 

0.39 

3.916 

709 . 1 

100.68 

0.40 


* Time taken assuming T - 0 for distance solidified = 25mm 


** y = location with respect to mid-plane 


Page 37 



TABLE 2 


RESULTS OF DECONVOLUTION USING THERMAL FIELD 2 



x ,mTn 

Temperature , °C 

y ,mm 

** 

y ,mm 

Time 

-3.916 

721.8 

26.15 

-1.24 

1 Hour 

-3.440 

720.5 

26.05 

-1.34 


-2.919 

718.5 

25 . 94 

-1.45 


-2.252 

716.0 

25.75 

-1.64 


0 

714.8 

25.74 

-1.76 


2.252 

714.5 

25.78 

-1.71 


2.919 

714.8 

25.88 

-1.62 


3.440 

715.7 

25.94 

-1.56 


3.916 

718.3 

26.11 

-1.39 


Time 

-3.916 

719.2 

32.87 

-1 ,34 

10 Hours 

-3.440 

717.6 

32.77 

-1.47 


-2.919 

716.2 

32.67 

-1.56 


-2.252 

714.0 

32.49 

-1.74 


0 

711.9 

32.43 

-1.91 


2.252 

712,2 

32.51 

-1.82 


2.919 

713.3 

32.65 

-1.69 


3.440 

713.5 

32.68 

-1.65 


3.916 

718.9 

32,84 

-1.50 


Time 

-3.916 

711.7 

101.06 

-1.68 

100 Hours 

-3.440 

710.2 

100.93 

-1.81 


-2.919 

709.0 

100,85 

-1.89 


-2.252 

707.9 

100.71 

-2.03 


0 

706.2 

100.55 

-2.19 


2.252 

707.9 

100.70 

-2.03 


2.919 

706.6 

100.74 

-2.00 


3.440 

708 .6 

100.86 

-1.89 


3.916 

709.2 

100.95 

-1.79 
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yet the interface is concave. For thermal field 2 the convex 
isotherms lead to a much more pronounced concave curvature in 
the interface. 

We also note from Tables 1 and 2 that the interface curvature 
changes slightly as solidification proceeds along the crystal. 

There is some scatter in the data due to the asymmetry of the 
compositional analysis and the drop in concentration occurring 
at the right hand side (positive x values) as seen in Figure 9 . 

The important point is that in both cases the interface curvature 

i 

is concave and is far less curved than what would be expected from 
the curvature of the 1 so concen t r a t ion lines. The interface does 
not follow the curvature of the isotherms as indicated by a con- 
cave interface consistent with convex isotherms in thermal field 
1 . This example indicates the importance of interface curvature 
on solute redistribution. A complete analysis using iteration 
between a thermal model and analyzed data is necessary to determine 
the actual degree of interface curvature encountered in HgCdTe grow- 
th experiments . 
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5.0 TWO DIMENSIONAL SEGREGATION 


5 . 1 BACKGROUND 

The results of crystal growth of HgCdTe shown in the previous 
chapters have illustrated some of the problems associated with 
the growth of this material. While longitudinal segregation 
would be expected as a natural result of directional solidification 
of a solid solution material, there is no obvious and simple ex- 
planation for the observed radial segregation. 

As has been discussed earlier, the problem can be considered 
In two parts. The first part is concerned with diffusion in 
2-Dimensions and will attempt to elucidate some of the work of 
other researchers and apply their results to HgCdTe. The second 
part is concerned with developing and demonstrating some of the 
equations and data for determining the effects of fluid flow 
and assessing the possible effects of "Shape Instability". 

5 . 2 TW O - DIMENSION A L DIFFUSION 

Lateral solution segregation has been the subject of recent papers 
2, 3 and it was felt that it should be carefully appraised 

to determine its applicability to HgCdTe. The interface compos- 
ition, and shape as determined by the deconvolution analysis de- 
scribed in the nrevious chapter can be incorporated in the model 
of Corieli et al (3) .They examined the unidirectional solidification 
behavior of a solid solution material with a curved liquid-solid 
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interface. The nomenclature, including figures, tables, symbols 
and equations is identical to those in their publication. 


Let V-*? velocity of solidification in +y m direction 

L=diameter of ampoule 

D=Diffusion Coefficient in liquid (D=0 in solid) 
x’^direction perpendicular to growth direction. 


Ampoule walls are at x f = t L/2 

Shape of the 1 lq u id / so 1 id interface is given by: 

y * = W 1 ( x ’ ) ' (5.1) 

Introduce the dimensionless parameters 

x = x T /L i y=y T /L, 6=VL/D and C-C’/C 1 ^ (5.2) 

The relevant equations can be written as: 

V 2 C+ P(3C/3y) = 0 (5.3) 

BC (k-JL) = ( 3C / 9y ) j - (ac/9x) t (3W/3x) (5.4) 

C ( x , y ■*■**>) = 1 (5.5) 


= 0 a t x = ± h • 


(5.6) 


C T is the concentration of solute in the liquid at the interface 


and C * oo is the bulk solute concentration. Subscript I refers 


to the interface. 


Coriell et al . assumed the interface shape to satisfy the 
equation W(x) = 6c.os (2irx), where cS is the amplitude across the 
diameter . This presents a problem in comparing our data to 
Coriell ’ s results in that the interface shapes we observe and 
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obtain from the deconvolution process are concave to the liquid, 
and the expression used by Coriell results in a convex interface 
for positive values of <$ . Due to the boundary condition at the 
ampoule wall, these two interface shapes may lead to different 
values in the change in solute concentration across the interface 
for the same deviation from planarity of the interface. It is 
instructive, however, to compare results for the same absolute 
values of deviations in composition and deviations from planarity 
for our concave interfaces and Coriel] Vs result for a convex 
interface. Further work should be done to obtain results for 
the same sense of interface curvature. With this restriction 
in mind, we may compare results using Coriell f s tabulated re- 
sults in their paper in Table 1 and their Figures 2-8. 

The crystal described in section 4 was chosen for comparison. 

The pull rate was 0.07 6 2cm/hr with a diameter of 0-8cm. Assum- 

-5 2 

ing an effective diffusion constant of 7x10 cm /sec we have 


and 

over 
ed ge 


^ / 4 


0. 0762 x .8 
7x 10" 5 x3 600x4(1 


= 0.019, 


(5.7) 


we take C 1 = .2 and k = 4 . 

00 

Coriell et alVs results correlate 4 , the interface amplitude 
a diameter, to the concentration at the interface center and 


(C 0 . (0) and C (- 2 )) respectively. 

v) 1 U I. 
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Table 1 in Cor iell * s paper includes values calculated for 
k = 0 . 4 , k=2 , and k= 1 0 . These were plotted and values for k = 4 

interpolated as follows; 


k 

3/4tt 

4 

C SI<°> 


4 

.01 

0 . 1 

1.037 

0 . 965 

4 

, 1 

0.1 

1 . 180 

0.745 


Note that the solute concentration at the center of the 
interface, 0^(0) , is higher than that at the edge. / This is 
due to the assumption of a convex interface. A concave inter- 
face will yield solute concentrations that are higher at the 
edges than at the center. As discussed before we will compare 
the magnitude of the concentration change to the magnitude of 
the deviation from planarity, although this comparison between 
concave and convex interface shapes may not be valid for larger 
deviations from planarity. 

Further interpolation was used to estimate values for 3/4= 

0.019 . 


k 

4 


3 / 4 tt <5 

0.019 . 1 


c si<°) c SI (%) 

1.05 0.94 


AC SI 

0.11 


The 

profiles 


crystal considered in section 4 with the concentration 
given in Figure 9 gives the following normalized con- 
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centration differences for an interface shape deconvoluted using 


thermal 

field 1 (see 

Table 1 

in section 4) 

* 





C SI <0) 

= 1 







c si ( -« 

= 1.14 

6 = 

. 19 

8x2 

= 0.012 




c sl ( + h) 

= 1.075 

6 = 

.08 

8x2 

= .005 



We 

see that even though 

the measured 

interface f rom 

section 

4 is asymmetric, we 

find AC 

values between . 
b X 

075 and . 

14 

cor res- 

pond ing 

to 6 values 

be tween 

.005 and .012. 

Based on solute 

t r an s- 

port by 

diffusion, Coriells 

results indicate that a AC - - 

O J- 

11, 

intermediate to our 

results. 

corresponds to an 

interface 

amplitude 

of 6 =. 

1 , which is 

an order 

of magnitude 

highe 

r than that 

ob- 


tained from the interface deconvoluted in section 4. Although 
this comparison is quite tenuous due to the number of assumptions 

made, we may make the following observations. 

Let us consider the maximum possible deviation from planarity 

for the interface. We argue as follows. For k>l, solute (CdTe) 

is preferentially incorporated from a bulk liquid of concentration 

C at steady state. For a concave interface, the solute concentr- 
o 

at ion will be highest at the edges and lowest at the center for 
either diffusion or shape instability. The minimum solute concen- 
tration possible in the liquid is obviously zero and the maximum 
possible liquid concentration is C • Thus the maximum solute diff- 
erential across the interface by any mechanism is (C o “0)“C o , and 
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thus the solid concentration difference is k C at a maximum, 

o 

HgCdTe liquid of concentration C q freezes at 805 °C , and pure HgTe 
(C = 0) feeezes at 668 °C , a temperature difference of 1 3 7 ° C . Under 
this extreme example, the center of the interface lies at a point 
1 3 7 °C cooler than the edge. If all isotherms were planar, this would 
correspond to a center to edge deviation from planarity of. 68cm in 
a 2 00 °C / cm gradient. Note that this is not necessarily the maximum 
interface curvature since we are stating that the interface edge lies 
1 3 7 ° C off of the 6 68 °C isotherm, curvature of the isotherms may re- 
sult in interface displacements more or less than given by planar iso- 
t he rms . 

i 

The concentration data for the crystal analyzed in our table 1, 

section 4, gave values for the solid concentration at the center of 

approximately C (0)=.2 and at the edge C (-%)=. 228, corresponding to 

s s 

a solidus temperature difference of about 5°C. This indicates that 
the maximum deviation of the interface from an isotherm for a gradient 
of 2 00 ° C / cm is 0.25 mm ( 5 ° C / 2 00 ° C / cm) . The observed interface curvature 
may be more or less than this due to curvature of the isotherms. Con- 
sider a constant gradient of 200 ° C throughout the solidification re- 
gion. If the interface lies in a region of convex isotherms (as in 
thermal field 1, section 4) but remains concave, a temperature diff- 
erence between the sample center and edge of 5 ° C will result in an 
interface curvature of less than 0.25mm. On the other hand, if the 
interface is convex, its deviation from planarity will be 0.25mm less 
than that of the isotherms. In a region of concave isotherms the 
deviation from planarity will be 0.25mm greater than the isotherms 
and the interface will be concave. 
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Returning to Coriell Vs results we find that to obtain a 

concentration difference of the order of AC /C =.11, we 

u l o 

must have a <$ = . 1 , which for an 8mm diameter ampoule used in 
the crystal discussed results in an interface deviation of 
1.6mm. Although such deviations are possible for this crystal, 
the interface would have to have been located in the cold zone 
with concave isotherms. The isotherms would have to be curved 
such that the center of the isotherm lay 1 . 6mm- . 2 5mm% 1 . 4mm 
below the same isotherm at the sample edge. Once again, this 

f 

is a possibility, thus more accurate measurements and improved 
thermal models must be used to ascertain whether lateral diff- 
usion alone can account for the radial solute distribution. 

We may state that for the reasonable thermal fields assumed 
here, the interface curvature is far too little to yield the 
observed solute distribution due to diffusion alone. 

Since solution of the two dimensional diffusion problem is 
of major importance in determining the cause of radial solute 
profiles, S EMTEC has developed its own two. dimensional diffusion 
model to calculate interface concentration for an assumed inter- 
face shape. We see to solve the same problem as Coriell e t al 
( 2 ) , namely 


V ? c + (V/D) ~ = 0 


(5.8) 


where C is the liquid solute concentration, V the translation 
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rate, and D the diffusion constant in the liquid. The coordinates 
x , y are centered at the sample mid-line with y positive into the 
liquid and the sample walls are at x = ±a. The boundary conditions 
are 


V ( k- 1 ) C j--D .[ ( 3.C / 3 y ) -(3C/3x) T (:3W/3x) ] , (5.9) 

C (y->oo) = C o , (5.10) 


and 9C/3x = 0 at x = ±a. (5.11) 

Here I denotes evaluation at the solid- liquid interface and the 
interface shape is assumed to be y=W (x) . 

Since this is a time independent problem, we may write 
equation 5.9 in a central finite difference form appropriate 
for numerical calculations. Let C. denote a point in the liquid 
not on the interface. Define 


h l 

ii 

X i-1 


h 2 

X i+ 1 

- X . 
X 


g l 

= y 3 ‘ 

y j-l 

(5. 12) 


g 2 = y j+r y j 

as the internodal distances from node (i,j) to its neighboring 
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nodes. After much algebra we obtain the finite difference form 
of equation 5.8 as 


C . . 

iJ 




(5.13) 


with 

P = 4/H 2 + 4 / G 2 - 2(h 2 -h 1 )(h 2 - h 2 ) /(AH 2 ) 
-2(g 2 -g 1 )(g 2 - g 2 > / ( G 2 B) + V ( g 2 - g 2 ) / ( B D ) 

Q=2/H 2 - 2h 2 (h 2 -h 1 )/ (AH 2 ) 

R=2 /H 2 + 2h 2 (h 2 -h 1 )/(AH 2 ) 

5 = 2/ G 2 - 2 g 2 (g 2 -g 1 )/(BG 2 )+Vg 2 (BD) 

T=2/G 2 + 2g 2 (g 2 -g 1 )/(BG 2 ) - Vg 2 /(BD) 


with 

2 2 2 
H = h* + h 2 



-f 


g 


2 

2 


(5.14) 

(5.15) 

(5.16) 

(5.17) 

(5.18) 

(5.19) 

( 5 . 20 ) 


A = h 2 h 2 + h 2 h ^ (5.21 ) 

B = 8ig 2 +g 2 8 l (5.22) 
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This rather complicated form allows variable node spacing 
in all directions to allow for chosing a node structure that 
gives the most accurate results. 

The boundary condition at the interface creates a problem 
in that solute concentrations are only defined for points on 
and in the liquid side of the interface. A central difference 
formalism is inappropriate here and we must write equations 
specific to either a concave or convex interface. Choosing a 
concave interface, we must use the backward difference form 

i 

for (3c/3x) , and the forward difference form for (3C/9y)^. 

Label the interface points by the index m and consider 

points along the interface with the same x values as chosen 

for nodes in the liquid. If the liquid, nodal points are at 

x.,v. then the interface points are at where x .“X 4 and 

x 3 i j 11 

Y . =W (x . ) . Def ine 
J i 


Cm = concentration of m-th point on interface, 

coordinates (x , W (x ) ) . 

m m 

CA = concentration of point in liquid above point 


interface. Coordinate 


y j+i) 



W(K m ) 



CL = concentration of point in liquid at (x^ 

W(x )), i.e. to the left of the interface 
m 

point for a concave interface. This value 


must be interpolated from concentration at 
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the nodes y. 11 > W (x )> v . . 

J + l- ™ - .3 

We may write 

ri C 

) I = (CA-Cm)/(y.-Y.), (5.23) 

dC/dx) l « (Cm-CL)/(x m -x rn _ 1 ) , (5.24) 

This gives the finite difference form of the boundary condition 
in equation 5.9 as 

C m = u [cA/ ( y j - Y j ) + c LOW/8x ) / ( x m -x in _ 1) ] (5.25) 

with u = V(k-l)/D+l/y.-Y.) + OW/3x)/(x m -X m _ 1 ). (5.26) 

Equations 5.13 and 5.26 form the basis for a finite difference 
method of solution of the two dimensional diffusion problem. Pro- 
gramming of these equations together with treatment of the sample 
edge and sample end boundary conditions has been underway. The 
total program is to date incomplete due to changes in the method 
of treating the interface nodes. Completion of the program for 
two dimensional diffusion will be carried out early in the next 
years effort. 
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5.3 SHAPE INSTABILITY 


The phenomenon of shape instability has been qualita- 
tively described in a previous section of this report. In 
this section an attempt is made to formulate this instabili- 
ty. No attempt is made to describe the most general problem, 
the interest being more in demonstrating the general charact- 
eristics and more importantly providing a theoretical basis 
for experimentation. Thus several simplifying assumptions 
are made while retaining the general character of the problem. 
As a first step we assume the following. 

A. The fluid motion is incompressible, that is 
the fluid density is not a function of the 
pressure (P) . 

B. The fluid density is a linear function of 
temperature (T) and concentration (C) . 

C. Heat transport by fluid motion is ignored. 

D. The dynamic viscosity v , the thermal con- 
ductivity and the specific heats are con- 

stant. 

These assumptions define the so called simplified 
Oberbeck-Boussinesq problem. If in addition the stress 
in the momentum balance is ignored, we may write the 
following equations . 
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3T 

at 


(5 . 27A) 


+ u • VT= k t V 2 T, 


ae 

at 


+ 


U • 


VC = D V 2 G , 


(5 . 2 7 B ) 


4rz + ( u • V ) u = - l/p n vP+v V 2 u + g (p/p -1), (5.270 

o c. o o 


V • u = 0 


(5.27D) 


with the constitutive relations 


P= P Q [ l-a(T-T o ) - 3(C-C q ) ] , 


(5 . 27E) 


C (}->«>) = C q , 


( 5 . 2 7 F ) 


c si" k c i 


(5.27G) 


T =T +m , C T +m 0 C_/ 
e m II 2 I 


(5..27H) 


Here k ^ is the thermal dif f usivity , u the fluid 
velocity, v the kinematic viscosity, F a pressure, g 
the gravity vector, and C the solute (CdTe) concentration 
in the liquid. Equation (5.27E) relates the fluid density 
p to the temperature and concentration and equation (5.27H) 
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approximates the interface temperature as a function of 
composition C ^ . The coordinate system (x,y,z) is an. in- 
ertial frame fixed to the ampoule tip with the y-axis 
along the sample axis positive into the liquid. The 
boundary conditions for the thermal field T are usually 
specified in a laboratory (furnace fixed) frame which 
moves with some pull rate R with respect to the ampoule. 

We now transform these equations into a coordinate 
system fixed to the center of the solid-liquid interface 
which moves with a velocity v . Let the coordinates be 

with n positive into the liquid. Denote the trans- 
formed fluid velocity as a ( C , n > O > and the solute concen- 
tration as y (£ , n , C ) = C(x,y 5 z). 

The transformed equations become 

-It- + 0 -?Y = DV 2 y , ( 5 . 2 8 A ) 

o t 

-fr “ v +(a*V)o = -1/p VP " +vV 2 a +g (p /p -1 ) , (5 . 28B) 
d L o O o 


V *0 = 0, (5.28C) 

P=P 0 [l-a(T-T o ) - 8(y-Y q ) ], (5.28D) 

Y (n^ 03 ) = Y 0 » (5.28E) 

Y SI “ k Y j > (5-28F) 
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(5.28G) 


T 

e 


T + 
m 


Y-i + 



Note that the only change is that we pick up an additional 
acceleration term in equation ( 5 . 2 8B ) . We have temporarily 
ignored the thermal field and will express it in the original 
coordinate system. The boundary condition at the solid- 
liquid interface is 


V ‘ n ( Y S r Yj) = DVy n , 


(5 . 28H) 


where v is the actual growth velocity and n the unit normal 
to the interface. We have not assumed a planar interface, but 
rather have moved into a frame attached to the center of the 
interface (assumed symmetric) whose velocity is v. Thus v-v 
at the center of the interface but not necessarily at other 
points along the interface. 

Equations (5.28A) through (5.28H), together with the thermal 
field equation provide the general basis for the study of shape 
instability. Initially we will decouple the thermal field by 
assuming a linear thermal field with flat isotherms in order to 
demonstrate the essential behavior of these equations. Later in 
the analysis we may reintroduce the thermal field equation or 
arrive at the thermal field by using a thermal model and inter- 
face deconvolution. 
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The first step in investigating the nature of shape 


instability is to assume that there exists a steady state 
3 o 3 y 

such that xf and tt vanish. The existence of such a 
3 t d t 

steady state has not been proven experimentally, however 
there is no reason to assume a priori that the features of 
shape instability would violate such an assumption. Further 
simplification of the equations can be made by making assump- 
tions about the nature of the fluid flow due to shape in- 
stability. We must exercise caution to avoid ’’throwing away” 


terms that are essential to the overall process of shape in- 
stability. The reduced set of equations can then be used to assess 
the influence of gravity and thermal field on interface shape. 

The full set of equations under steady state conditions 
may be solved numerically given the proper boundary conditions. 
Note that we are interested in a solution of the equations in- 
volving macroscopic deviations of the interface from planarity 
under conditions that would normally lead to a thermo so lu tally 
stable configuration, i . e . , hot end up an d increas ing density 
downward with respect to gravity. This situation and the re- 
sulting interface shape instability have not been addressed 
by other workers investigating convective instabilities or 
interfacial instabilities. Conclusive theoretical evidence 
for shape instability will be a prime justification for solid 
solution crystal growth in a micro gravity environment. 


Page 55 



We may derive a generalized relation for the mass balance 
along the interface. We assume that a steady-state exist such 
that and ~ vanish. Furthermore we assume that in the 

u t d t 

steady-state the interface shape remains constant and thus the 


growth velocity v - v^rj and that the ampoule is sufficiently 

3 y 

long such that there exist an Hoo so that y*y and there 

o 3 n 

is no fluid motion due to convection as . Consider the 

following relations from the set of equations described earlier. 


(5.29) 


Z'Vy= DV ? y 


and 


v • n (k- 1 ) y j = D 7 y 


A 

* n 


If the interface is denoted by 


(5.30) 


(5.31) 


■r. = W ( £ ) (5.32) 


then the unit normal to the interface is given by 

* - B(^-(3W/3f.) £) (5.33) 

where 

B = ( l + OW/94) 2 )~ 1 ' 2 (5.34) 
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writing 


A A 

a = o £ +o 2 n (5.35) 

then equation (5.30) becomes 

Df 9 2 y/ 3C 2 + 3 2 y / 3 n 2 ] - o l (3y/ 3£) +o 2 (3y/3n) (5.36) 

This may be written in a different form as 

3/3£(3y/3£-o 1 y/D) + f- ( 9 Y / 3 q-o 2 y/D.) 

+ (y/D) ( 3 a j / 3£+9a 2 /3n) - 0. (5.3 7) 

Note that the boundary condition in equation (5.29) gives 

3 a | / 3 £ +3 o 2 / 3 n = 0 , (5.38) 

thus equation (5.37) reduces to the form 

3 N/ 35 - 3M/3n= 0 (5.39) 

with N = 3y / 3 £ -a j^y (5.40) 

M = -(3y/3n ~ CT 2 Y ) • ,(5.41) 

Now Gauss’ theorem is of the form 
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(5.42) 


//(3N/H- 3M/3n) d n d £ = fMd.£ + Ndn 
R C 

for a region R bounded by a closed path C. Integretion of eqn. 
(5.39) and application of Gauss' theorem gives 

f - ( 3 y / 3 h - a £ Y / D ) d £ + ( 3 y / 3 £ -o ^ y /P ) dp =0 (5.43) 

C 

Consider the sketch below representing a sample, interface 

' f " 

and the liquid region bounded by an ampoule at n=±a . 



We wish to evaluate equation (5.43) along the closed path 
A-B-C-D-A. Write the integral as a series of line integrals 


B C D A 
f = / ■+ / + / + / = 0 . 

C A B C D 

We note that along A-B we have the following 


;=a , dr=0 


_3T 
3 5 


= 0 , a - 0 


and thus 


B 

/= 0 . (5.44) 

A 


In a similar manner we note that 


D 

/ =0. (5.45) 


Evaluation of the path integral along B-C we note that 


n=noo > dn = 0 ,a 1 = 0,a 2 =-v o 


giving 


C 


/ 


B 


- a 

/ -( 3 y / 3 n ) v^*Y /D ) d 5 

a 


h= n 


oo 


(5.46) 
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a a 

/ = / ( ( Y / D ) ( a 2 “ ° i 3 W / 3 O + ( 3 Y / 3 £ ) ( 3 W / 9 £ ) - 3 y / 3 q ) d £ 

D -a 

(5 .48) 

where all quantities are evaluated along the interface. 

Combining all the path integrals gives 

a 

f [ y ( o 2 ~ o i 3 W / 3 r ) +D ( ( 3 y / 3 f, ) ( 3 W / 3 £; ) - 3 y / 9 q ) ] df^-2 av Q Y 0 

- a 

(5.49) 

Now we use the boundary equation at the interface, equation (5.31) 
t o write 

(k-.l ) v o y x = f) [ ( 3 Y / 3 C ) ( 3 W / 3 f ) - 3 y / 3 n ] (5.50) 
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Equation (5.49) then becomes 


where we have taken v « v H . 

o 

a 

/ (° 2 ” 0 1 ( 3W/.3.5) - (k~.l ) v 0 ) I d 5 = -2a o y o (5.51) 

- a 

This is the generalized solute balance equation where all 

quantities are evaluated along the interface. It relates the 

net solute balance due to segregation, motion of the interface, and 

fluid flow to the solute in the bulk liquid y . The equation 

o 

(5.51) reduces to a more familiar form if we assume /that the only 
fluid motion is that due to advancement of the interface (remember 
that o is the fluid velocity measured in the interface fixed co- 
c ordinate system). For this case we have 

Cj^O, o 2 = -v q giving 

a 

/ k V = 2a ' Y o ’ (5.52) 

-a 

which is the familiar form of solute balance along an interface 
at steady state. 
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6.0 SUMMARY AND RECOMMENDATIONS 


In the preceding sections of this report we have address- 
ed the characteristics of solute redistribution in the HgCdTe 
system from both a theoretical and experimental point of view. 

Various segments of the work have concentrated on both one- 
dimensional diffusion and two dimensional aspects of solute 
redistribution. The overall goal of this effort is to identify 
and quantify those factors resulting in compositional in homogen- 
eities in directionally solidified HgCdTe. 

The one dimensional model developed in a previous year’s work 

i 

has been applied to several experimental cases and the sensitivity 
of the model results to variations in thermal gradient and ampoule 
pull rate examined. Cases in which the ampoule pull rate was 
changed or in which growth was interrupted for a specific time were 
modelled to assess the feasibility of using rate changes to identify 
the interface location and the resulting axial solute profile. Al- 
though the utility of the model is limited by its one dimensional 
nature, it has obvious applicability in selecting certain funda- 
mental processing parameters (ampoule length, pull rate, etc.) to 
obtain an axial profile with a sufficiently long steady state region. 

Although axial solute changes are an unavoidable result of 
diffusion controlled growth, the presence of radial solute changes 
may be eliminated under proper growth conditions. The nature 
and causes of these radial inhomogeneities is to date unknown. 

We have addressed two main mechanisms for radial solute varia- 
tions; lateral diffusion due to a curved interface and shape 


Page 62 



instability caused by gravity driven flow along a curved inter- 
face. These effects are not mutually exclusive, certainly 
diffusion may take place in the presence of fluid flow. In 
both cases the interaction of the thermal field and solute dis- 
tribution in effect determines the interface shape which in turn 
depends on the driving force for the solute distribution in 
the liquid. Knowledge of the thermal field and interface shape 
is fundamental to the application of either a diffusion model 
or a shape instability model to actual experimental data. 

The technique of interface deconvolution provides a method 

i 

of determining interface shape and composition from analysis 
of a grown crystal and knowledge of the thermal environment. 

Several examples of interface deconvolution have been given to 
illustrate the consequences of an interface with a varying solute 
distribution being neither an isoconcentrate nor a isotherm. 

As thermal models are augmented by experimental measurements 
and improved values of the thermophysical properties, the technique 
of interface deconvolution will play a large role in determination 
of the interf ace shape and composition. These data are vital to 
the understanding of the causes of compositional variations. 

The mechanism of lateral diffusion along a curved inter- 
face has been reviewed and comparison to results obtained by 
other workers has been made for a HgCdTe crystal. Although 
no definite col elusion can be made at this time, it appears 
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that the degree of interface curvature necessary to give the 
experimentally observed radial solute variations is far greater 
than the curvature that would be expected in normal processing 
conditions- More refined work in this area is needed. 

Finally, we have indicated appropriate equations and boun- 
dary conditions for the problem of shape instability. A gener- 
alized interface solute balance was derived in the presence of 
fluid flow. These equations and conservation relations form the 
basis for perturbation calculations to determine the extent of 

i 

shape instability. These calculations are to be carried out 
in the following year. 

In order to address all of these issues, certain exper- 
imental and theoretical techniques must be applied. Knowledge 
of the thermal field is paramount to a detailed interpretation 
of the solute profile . Although we may proceed theoretically 
using an assumed thermal field and investigate the qualitative 
behavior of the system, comparison of experiment to theory 
depends on a detailed knowledge of the actual thermal profile 
in the crystal. Accordingly, characterization of the growth 
furnaces should be a first priority. Since it is difficult, 
if not impossible, to instrument a HgCdTe ampoule, materials 
of similar thermal conductivity may be used to profile the 
furnace. At the same time, efforts to obtain the best values 
for the thermophysical properties of HgCdTe and their temperature 
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temperature dependence should be undertaken. The furnace pro- 
files can then be used to fine tune the existing thermal models 
which can then be applied to the HgCdTe growth system. The 
goal of this effort should be to obtain the best estimate of 
the thermal profile at steady state. 

Analysis ^of the resulting solute profile in the grown 
crystal is of obvious importance. Use of the SEM with an energy 
dispersive spectrometer provides the only method with adequate 
spacial and compositional resolution. Further work with WDS 
and improved standards will yield higher absolute accuracy in 

f 

the x values. Computer assisted analysis of the composition 
profiles will reduce the time for a complete composition profile. 
The analysis of the grown crystal coupled with the thermal field 
can be used to deconvolute the interface shape and composition 
as described earlier in this proposal. 

In order to compare the results of this interface decon- 
volution with observed values, we must be, able to delineate the 
solid-liquid interface at some point in the growth. Peltier 
pulsing of the interface is a very renote possibility due to 
ampoule instrumentation problems. More likely approaches would 
depend on a rapid quench, decanting, or perturbation of the solute 
profile by a rate change or shaking of the ampoule. Each of 
these techniques has distinct advantages and disadvantages. 
Experiments that compare interface shapes at the same point in 
growth under the same thermal conditions should be carried out 
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with the interface marked by a quench, a rate change, and a 
physical disturbance (vibration, etc,) • These interface shapes 
and compositions may then be compared to those obtained by inter- 
face deconvolution. 

The effect of ampoule design and early growth conditions 
must be quantified. Previous experiments have used a conical 
shaped tip on the ampoule to promote grain selection. Other 
techniques, such as a neck-down region should be assessed. 

The extent of thermal coupling plays a large role in determ- 
ing the initial interface shape. Both experimental and calcul- 
ation techniques should be used to determine at least the pos- 
ition in the furnace of the first to freeze point. 

All of these experiments have been directed towards 
quantifying the influence of growth parameters (temperature, 
pull rate, etc.) on the solute profile in the solid. Know- 
ledge of the interface shape and composition is important from 
a theoretical standpoint if we are to ascertain the role of 
shape instability. 

There are additional experiments that address directly 
the role of gravity in producing compositional inhomogen it ies . 

By a suitable mechanical arrangement, the growth furnace may be 
tipped at an angle with respect to the g-vector while keeping 
the ampoule at the center of the thermal field. Crystals grown 
under these conditions should show a distinct asymmetry in their 
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solute profile and interface shape if gravity driven fluid flows 
are important. Convective instabilities not related to those 
causing shape instability may be present in this case however. 

Another experiment represents the extreme limit of the one 
discussed above. By careful ampoule design, it may be possible 
to grow in the unstable mode (hot end down) and thus promote 
thermo-solu tal convection rather than suppress it. The value 
of this experiment should be considered in light of the certain 
presence of bulk convection which may obscure the effects of 
shape instability. 

Finally, we should seek to discover model systems for HgCdpe. 
These would be characterized by a material that shows a widely 
separated liquid. us and solidus, rejected component more dense 
than the bulk liquid, and hopefully transparent to visible light. 
The essence of shape instability could be best approached ex- 
perimentally with an ideal model system. 

For any set of experiments, the interplay between theory 
and experiment will follow the general sequence outlined below. 

1 . Characterize furnace using instrumented 
ampoule and model material . 

2. Apply thermal model to estimate thermal pro- 
file to be obtained in HgCdTe ampoule. 

3. Use one dimensional model to select growth 
rate and ampoule length giving the best com- 
promise between length of transient regions, 
thermal gradients, and ampoule length. We 
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desire an extended steady state region both 
compo s it ionally and thermally. 


4. Perform experiment 

5. Analyze the crystal composition profile using 
SEM/EDS. Reduce data into appropriate format. 

6. Use the interface deconvolution model and the 
thermal model iteratively to arrive at a self 
consistent interface shape and thermal field. 

■■f T 

7. The fluid flow equations in Section 5 can now 
be solved using the interface shape, composi- 
tion and thermal field to give the liquid 
composition and fluid velocities as a function 
of g. 

If the resulting fluid velocities and liquid composition 
are consistent with those postulated for causing shape instab- 
ility, we can conclude that shape instability is a dominant 
mechanism for creation of radial in homogeneities and that growth 
in a microgravity environment should eliminate these inhomogen- 
eities. 

In conclusion, we have addressed the key issues and ex- 
periments necessary to quantify the causes of solute redistrib- 
ution in HgCdTe . The results of these experiments and analysis 
will allow one to identify the extnet of influence of processing 
parameters on the solute profile. By combining experiment and 
theory we may then ascertain the degree to which gravity influences 
this solute redistribution. The ultimate test is then to perform 
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experiments under microgravity conditions in space. 
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IN DIRECTIONALLY SOLIDIFIED HgCdTe 
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ABSTRACT 


A one dimensional model for solute redistribution is developed 
that can be applied to solid solution alloys with widely separ- 
ated iiquidus and solidus phase diagrams. The model is applied 
to the HgCdTe system and accounts for the variation in interface 
temperature, segregation coefficient, and growth velocity with 
composition. The model considers a finite length ampoule, treat- 
ing both initial and final transient segregation. Agreement be- 
tween the model and experimental results is good, indicating that 
the axial solute redistribution in directionally solidified HgCdTe 
crystals is dominated by diffusion. Comparison of the model to 

experiment in the final transient region gives values for the 

-52 -52 

effective diffusion constant of 5x10 cm /s to 7x10 cm /s. 
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1 .0 INTRODUCTION 


One of the major problems in production of HgCdTe crystals 
suitable for large array infrared detectors is material homo- 
geneity. Axial and/or radial composition variations lead to 
variations in detector response over the area of the detector. 
Growth techniques that eliminate or minimize these variations 
are of primary importance in meeting the demands of high 
quality HgCdTe crystals. 

In this work we model the directional solidification grow- 
th method with emphasis on the HgCdTe system in an effort to 
understand those growth parameters influencing crystal uniform- 
ity. Cur approach is to treat the one-dimensional planar inter- 
face solidification problem with no convection, in a similar 
manner to Favier [1] and other more classical approaches [2,3] 
with certain important characteristics of the HgCdTe system 
taken into account. 

HgCdTe is representative of a number of alloy semiconductors 
which form solid solutions over their entire range of compositions. 
The behavior of the HgCdTe system is such that some of the assump- 
tions made in the studies mentioned above are violated. Figure 
1 is a representative pseudo-binary phase diagram for the H ; g ^ - 

Cd^Te system [ 4 ] . Consider an initial liquid with CdTe (solute) 
mole fraction x = . 2 . We then find the following characteristics 
of this system which are different from those found in dilute 
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alloy systems • 


(1) The interface temperature varies from almost 

800 ° C (T q ) initially to 705°C ('T ) at steady 

state and drops to 6 6 8 ° C over the final tran- 
sient region of solidification. 

(2) The 1 iquidus and solidus curves are widely 
separated and yield an equilibrium segregation 
coef f ic lent , k, that var ies significantly with 
liquid composition . 

(3) As a consequence of (2), the solid being formed 
from a liquid in equilibrium with it have widely 
different compositions. For example, the first- 
to-freeze solid from a z= . 20 mole fraction liquid 
has x = . 5 4 , and steady state results with a x~ 10 5 5 
liquid forming a solid of x=.20. 

The goal of this work is to develop a one dimensional 
model for directional solidification that takes into account 
the factors above which distinguish the HgCdTe and other solid 
solution systems from the dilute alloy or doped semiconductor 
systems treated by previous workers. This model will then be 
used to analyze the transient segregation experimentally ob- 
served in HgCdTe. 
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2.0 THEORY 


In this section we wish to formulate the most general 
one-dimensional solidification problem in keeping with the 
directional solidification ( Bridgman- St ockbarger) growth 
configuration and the KgCdTe solid solution system. Accord- 
ingly, we will make the following assumptions: 


(1) The thermal and concentration fields are one 
dimensional , 

F 

(2) No significant diffusion occurs in the solid. 

(3) Convective effects are negligible. 

(4) The thermal field is a specified function 
of the spacial variables and time and does 
not change as a result of the solidification 
process. 

(5) The sample movement rate, R, is constant. 

The interface velocity or growth rate, V, 
is not necessarily constant. 

We choose a coordinate system such that z is measured 
from a fixed laboratory frame (the furnace mid-plane) and 
x is measured in the sample from the ampoule tip ( first- to- 
freeze point). The variable n is the distance into the liquid 
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measured from the solid-liquid interface (see figure 2). 

Let us denote variables evaluated at the position of 
the interface by a star (*).. Choose t = 0 to be the start 

ft 

of solidification in the ampoule tip and z (0) the z~ co- 
ordinate of the ampoule tip at t =0 . The ampoule is being 
pulled downward at a rate R with R taken as a positive 
number. We may then relate the various co-ordinates by 

x = z-fF.t-z (0) 

ft ft 

n=x-Rt+z (C)-z (t) 

The total length of the ampoule is 9 . , and £ * = denotes 

the distance from the interface to the top end (last-to- 
freeze) portion of the ampoule. 

It is most convenient to express the thermal field in 
the fixed coordinate system, and we shall assume a thermal 
field of the form 

T ( z , t ) = T +Gz ( 3 ) 

o 

where G is the thermal gradient, assumed to be constant. 

Let L(C) denote the 1 iquidus temperature for material 
of composition C, and let C(n 5 t) denote the solute (CaTe) 


( 1 ) 

( 2 ) 
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mole fraction in the liquid ahead of the solid-liquid inter- 
face. The interface location at any point in time is in 
effect determined by the thermal field. The solute concen- 
tration at the interface is given by C(G,t) and the correspond- 
ing liquidus temperature is L (C (0 , t ) ) . We may then define 

k 

the interface location z such that 


T ( z * , t ) = L ( C ( 0 , t ) ) 


(4) 


The solute concentration in the liquid satisfies 


D 


3 : C 

3rT 


+ V(t) 


3C 
3 n 


3_C 
3 1 


(5) 


where D is the effective diffusion constant and V(t) is the 
growth velocity. The boundary conditions are 


C ( r, 0 ) = C 

o 


t = 0 


( 6 ) 


3 v ( n , t y _ g at ~ = i » f or a ij_ t 
3n 


(7) 


V ( t ) (k(C) - 1) C(0, t) = D t>0. 


( 8 ) 
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The boundary condition in eq, (7) is normally written 


for an infinite amnoule as, C ( n , t ) = C as for all 

o 

In this study we wish to solve for both the initial ana 
final solidification transients, hence we require that 
there be no source of solute at the ampoule end ( q = £ f ) 
giving the boundary conditions in eq . (7). 

In addition to the above equations, we shall assume 
that the liquid us portion of the phase diagram over the 
region of interest may be approximated by 

L ( C ) - T . +M A C+yC 2 
A A 

where T , and y are appropriate constants obtained by 

least-squares fitting. We may then write 

v(t) - a+|f 

* 

where z satisfies 

T ( z * , t ) = L ( C (0 , t ) ) . 


(9) 


( 10 ) 


(ID 



Using the assumed expressions for the thermal field and 
the liquidus curve we find 


z*(t) - t/<5(T A -T Q +K A C(0,t)+ Y C 2 (0, t) ) , 


giving 


dz 

dt 


(l/G) (K A +2yC(0,t)) — 


The numerical procedure used to solve equations (5 

to (8) is unique. Let the incremental time and special 

steps be At and An respectively. Assume at time t Q the 

values of C(n»t) have been determined for all t< t and 

— o 

all n . The goal is to then compute C(n,t T ) for all n 

* 

where t'~t + A t • We assume a trial value C for the con 
o 

centration at the interface at time t f , say 
C(0, t \)=c*. 


We then compute 

qC ( 0 , t T ) 
9 t 


C' A - C(0,t ) 

i 1 q^L 

At 


( 12 ) 


(13) 


(14) 


(15) 


F a « e ) 9 



giving the interface velocity 


V(t')=R+|f | t = t , =R+(1 /G)(K a +2 Y C*> • (16) 

We may then compute the first spacial derivative via the 
boundary condition (eq. 8) as 


3 C ( 0 , t ' ) 

3 n 


( 1 / D ) V ( t ’ ) ( k ( C * ) -1) C ( Q , t ' ) 


(17) 


where k(C ) it determined by the phase diagram and is in 
practice obtained by interpolation in a table of values for 
k versus C. The second derivative is found via the diffusion 
equation (5) giving 


3 2 C (0 , t ' ) . . , 3 C ( 0 , t 1 ) 

aTT 2 u / d ; { 3 1 


- V ( t * ) (18) 

9 n 


We then may obtain our desired results by expansion 


C (An , t 1 ) =C (0 , t * ) ~ An c - Q - ’ 2 — • 


An 2 (19) 


for the value of the concentration in the liquid a distance 
drj away from the interface. 
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We now propogate C(n,t) out in until the end of the 
ampoule is reached ( n = £ T ) and test to see if the boundary 
condition 


dC 

3 'n 


( n , t * ) =0 at n = £ f 


( 20 ) 


is met. The solution may be propagated to n 1 = 2 An , the next 
spacial step, by first calculating 


3 C ( A n > t 1 ) 
at 


, t 1 ) “C ( An , t _ ) 

— ■ • — — o~ 

At 


(2 1 ) 


We then use a Taylor series expansion and the diffusion 
equation to obtain 


.3C(Ah,t ’) 

a n 


3 C ( 0 , t ' ) ' 3 2 C ( 0 , t ' ) 

a n or" 


An * (22) 


and 


3 2 c ( & n , t r ) n ,3C (An , t ' ) 3C(An,t 'K 

" 3n 2 ’ (x/D)( 3t V(t ' 3n ' 

( 23 ) 
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second expansion leads to the desired results 


c(2An,t')=c(An,t’)+ 


3C(An,f) 
3 n 


, 3 2 C ( A n , t * ) An" 

3.n ‘ 2 


(24) 


The same process is used for the next spacial step 3An, etc. 

until n = £'» the end of the ampoule at this time step. A test 

is made to see if the boundary condition eq . (20) is satisfied. 

* 

If it is not, a new trial value of C is selected and the entire 

' i ' 

process repeated until the boundary condition is satisfied. 

The initial value boundary condition 

c (n ,0) -C Q (25) 

allows one to take t Q =C and proceed in a straightforward manner 

to calculate C(n, t) for any t. At each time step, after the value 
* 

of C(0,t)=C has been found that satisfies the boundary condition 
in eq. (20), we may find the solute concentration in the solid as 


C (>:*) - k (C ( 0 , t ) ) C(0,t) 
s 


(26) 


where we note that 2 (t), and thus x ( t ) , are determined by 

the thermal field and i iquidus curve of the phase diagram. The 
remaining amount of liquid is given by 


baye g2 



(27) 


£ 1 


* 

£-x 


The time is incremented by At, and the entire process re- 

* * * 

peated to yield new values of x (t), z (t), V ( t ) , and C (x ). 

s 
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3.0 RESULTS 


In this section we present results from the model and 

compare these results to experimental solute profiles obtained 

from directionally solidified HgCdTe crystals. In figure 3 

we present a typical solid concentration as a function of 

position taken for an initial liquid composition of x= . 2 . 

Theoretical curves of solid composition are plotted for two 

_ 5 2 

different values of the diffusion coefficient D-5x 1 0 cm /sec 
-5 2 

and D=7xl0 cm /sec. We note the initial transient region 
marked by a decreasing solute concentration, a steady state 
region extending some 105mm, and a final transient in which 
the solute concentration tends to zero as the end of the 
ampoule is reached. A solute distribution of this nature is 
typical for a segregation coefficient k > l , and in fact k varies 
from 2.7 initially to 3.6 in the steady state region. At 
this translation rate, the transient regions are rather short, 
and the differences between the curves with different diffusion 
constants are most marked in the final transient region. 

In figure 4 we plot both the growth velocity and inter- 
face position relative to the furnace as a function of distance 

along the crystal. The other input data is the same as in 

-5 2 

figure 3 and we chose 0=5x10 cm /sec. The thermal field is 
assumed to be such that the steady state solidification 
temperature (705°C) occurs at the furnace midplane (z-0). As 
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a consequence of the phase diagram, the first to freeze from 
a x=,20 liquid occurs at a higher temperature (%800 ° G) and 
consequently is located higher in the furnace. As solute 
is preferentially incorporated, a diffusion layer is estab- 
lished and the liquid composition in front of the interface 
approaches the steady state value of x= . 055 which solidifies 
at a temperature of 705°C. During the approach to steady 
state, the interface position changes, moving toward the 
furnace midplane. Thus the interface moves along with the 
ampoule toward the cooler portion of the furnace, although 

f 

not at the ampoule pull rate R . At steady state, we have 
a fixed solidification temperature and hence a fixed inter- 
face position in the furnace. As the end of the ampoule 
is approached, there is no longer enough solute remaining 
in the liquid to replenish the diffusion layer, resulting in 
a decreased solute composition at the interface. This, in 
turn, results in a still lower solidification temperature and 
consequently the interface position moves toward the cold end 
of the furnace. 

The growth velocity, or the velocity of the interface 
relative to the ampoule tip is also shown on figure 4. As 
discussed above, the interface moves down toward the furnace 
midplane during the initial transient region, resulting in a 
growth velocity that is initially zero and increases until the 
steady state region is reached. At steady state the interface 
is fixed with regard to the furnace and consequently the growth 
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rate and ampoule pull rate are equal. The final transient 
region presents some rather unexpected behavior. As the solute 
is exhausted in the remaining liquid during the final transient, 
the solidification temperature decreases and the interface again 
moves along with the ampoule toward the cold zone, resulting 
initially in a growth rate less than the translation rate. At 
the very end of the ampoule, both the solute concentration and 
the concentration gradient in the liquid approach zero due to 
the boundary condition at the ampoule end. As a consequence, 
the time rate of change of solute concentration at the interface 
also approaches zero, resulting in the velocity increasing back 
to the ampoule pull rate at the ampoule end. The numerical 
scheme used here may over-estimate the length of this region 
due to a finite node size, but the effect is present, though of 
little technological importance. 

We have indicated some features of the model as it applies 
to the HgCdTe system. Figure 5 shows a plot of experimentally 
measured solute c on cen t r a t ion s and calculated solute profiles 
from the model. The data is taken from a HgCdTe crystal with 
C o =.20 and an estimated thermal gradient of 200°C/cm. In this 
instance the translation rate was .762mm/hr. The experimental 
data points were determined by using an SEK with an energy dis- 
persive spectrometer and HgCdTe standards to permit peak de- 
convolution [5] . Since radial concentration variations occur in 



this system, the axial composition profile is an area-weighted 
average of the radial compositions arms s oar h s Ti re at- a 
particular axial position. The results of such a procedure 
would be equivalent to those obtained by density measurements . 

There are several features to note in figure 5. Agreement 
between the model and data in the initial transient region is 
not very good. This may be due to several effects. Initial 
solidification occurs in a conical ampoule tip, conditions which 
do not approximate those assumed here. Secondly, although not 
shown on this graph, the experimental data shows an -initial 
increase in solute concentration in the tip of the crystal be- 

- r : i * :i j A qr 

fore it begins to decrease. This is thought to be due to under- 

- - .. ■ f. 

cooling which certainly must occur to some extent and can lead 

r 

to nucleation away from the tip. The net result of this is 
to cause the initial transient region to occur further along the 
crystal then anticinated. resulting in disagreement with the 

calculate a curves. 

1 - . | 

■ ■ - “ .i - w .... - i J •-} V ■'■■■> r 

In the steady state region, agreement between experiment 
and theory is quite good as would be expected. The data lies 
below the calculated curve as the final transient region is 
approached. Although the experimental analysis has a resolution 
od Az'vi .002 , the availability of absolute standards my limit 
the absolute accuracy of these measurements. In addition, in 

the model computational time constraints result in chosing a spa- 

# 7 

cial and time step size sufficiently coarse as to yield accumulated 
truncation errors in the later stages, giving a steady state 
value of X-. 205 instead of x= .2 00 . Taking both of these effects 

UU 
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inee g©n§id§rg6i©n, eh@ agriiininE i§ gg©d fgr a m©i§l with n§ 
feadfea§ic ffgffl the §xp§r im§ntal data. 

The final eranaiant ragian is best axaminad eg an expanded 
seals, figure I is an enlarged pgrtlan gf the final transient 
regien. In this ease the experimental data have been shifted 
upward and the calculated curve downward to yield a value of 
in the steady state. In addition, the ealeulated eurve 
has been shifted toward the origin very slightly to account for 
errors in position measurement. The same corrections were 

i 

applied to all data points shown, hence we are in a sense ad- 
justing the data to conform to the calculated curve but have not 
changed the shape of the data curve. It is apparent that the 

shape of the data curve most closely fits a value of the effective 

« 5 2 

diffusion constant D»7xl0 cm /see. We note that the shape of the 

-5 2 

ealeulated eurve for D»5xlQ cm /see is not eonsis tent with the 
data even if we translate the position of the data points. In 
other crystals grown under different thermal fields and trans- 
lation rates, we have found that the experimental final transient 

- 5 2 

eurve is always fit using values for D between 5x10 cm /see 
and 7xl0 m ^em^/see [6J. 

4, CONCLUSIONS 

We have presented a one dimensional model for calculation 
of axial solute distribution in the HgCdTe solid solution system. 
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The model accounts for the change in interface temperature, 
segregation coefficient, and growth velocity due to a widely 
separated liquidus and solidus phase diagram. The model also 
accounts for the effects of a finite length ampoule, allowing 
the final transient segregation to be calculated. 

As a consequence of the phase diagram, we see that the 
interface velocity and position vary over the initial transient, 
remain constant in the steady state, and vary again over the 
final transient. The implications of this are that obtaining 

flat solidification isotherms over the entire region of growth 

i 

in a real furnace system can be difficult if we must also main- 
tain a high thermal gradient. 

By comparing experimental results to calculated curves 

in the final transient region, we conclude that the effective 

- 5 2 

diffusion constant in the HgCdTe system lies between 5x10 cm /sec 
- 5 2 

and 7x10 cm /sec for an initial liquid composition of C q = . 20 . 

Agreement between experiment and theory is good indicating 
that diffusion dominates the axial segregation in the HgCdTe 
system. In a subsequent paper, we will treat the two dimensional 
problem and interface stability. 
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FIGURE CAPTIONS 


FIGURE i 


FIGURE 2 


FIGURE 3 


FIGURE 4 


FIGURE 5 


FIGURE 6 


Pseudo-binary phase diagram of Hg ^ ^Cd^Te . 

Co-ordinate system used in model. The ampoule 
* 

tip is at z (o) at t = o. 

Calculated solid axial solute concentration 

versus distance along the crystal. Curves for 

-5 ? -52 

D=5x 10 cm"/s and 7x10 era / s are shown . 

- Interface velocity and location relative to 
the furnace as a function of distance along 
the crystal. The scale is broken to eliminate 
the constant steady state region. Input par- 
ameters to the model are the same as in figure 3 . 

Comparison of experiment and model results for 

a directionally solidified crystal. The initial 

liquid composition was z = . 2 0 . Calculated curves 

-5 2 ' _ 5 2 

are shown for D=5x 10 cm /s and D-7xl0 cm /s. 

Enlargement of the final transient region of 
figure 5. The data and calculated curve have 
been shifted to yield the best agreement (see 

text) , The upper curve is calculated with 

-5 2 -5.2 

D = 5xl0 cm / s and the lower curve with D = 7x 10 cp / s . 
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INDIRECT DETERMINATION OF THE INTERFACE 


IN 

SOLID SOLUTION CRYSTALS 


Donald C. Gillies and James Creed Clayton 
Semtec Incorporated, Huntsville, Alabama 35805 


ABSTRACT 


A technique is described by which the location, shape and 
solid composition of the solid-liquid interface in a solid- 
solution crystal is determined by indirect means. The method 
utilizes experimentally measured composition data from the 
solidified crystal together with thermal profiles calculated 
from a thermal model to iteratively determine interface par- 
ameters consistent with both thermal boundary conditions and 
the measured solute distribution. This technique allows for 
an interface that may be neither an isoconcentrate nor an iso- 
therm. Typical interface shapes in HgCdTe crystals derived by 
this technique are given as examples. 
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1 . 0 INTRODUCTION 


Solid solutions are presently being utilized in several 
semiconductor applications with promise for wider applicability 
in the future. Mercury cadmium tellur ide (HgCdTe) , for example, 
shows superior performance as a radiation detector in the short 
to medium infrared region with potential for extension to long- 
er wavelengths. The directional solidification of solid solutions 
has, however, met with limited success. Apart from the chemical 
and structural problems of a particular system, directionally 
solidified solid solutions show pronounced compositional in- 
homogeneities. These inhomogeneit ies limit the materials yield 
for single devices and the performance of large area arrays. 

These compositional inhomogeneit ies, present to some degree in 
most solid solutions, depend to a large degree on the shape of - 
the solidifying interface. The determination by experimental 
and theoretical means of the interface shape during the undirect- 
al solidification of a solid solution system is the subject of 
this paper. 

In the HgCdTe system the high degree of preferential in- 
corporation of CdTe leads to severe inhomogeneities . Proper 
choice of ampoule length and crystal pull rate will result in 
a steady state growth region with a constant axial solute com- 
position, however, compositional analysis of crystals indicate 
that there still exists a high degree of radial segregation. 
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The cause of this radial segregation is largely unknown. 

Possible mechanisms include lateral diffusion across a curved 
interface (1,2) convective effects and thermal effects (3*4). 
Another postulated mechanism in HgCdTe is a result of the 
heavier rejected component (HgTe) flowing along a curved inter- 
face leading to higher concentrations of HgTe at the lowest 
points in the interface. If the interface were concave, the 
increased HgTe concentration at the center of the interface 
would lead to a lowering of the solidus temperature at the 
center relative to the edges, further increasing the curvature 
of the interface. This global interface instability has been 
termed shape instability by Davidson ( 5) and has been observed 
in other systems in which the rejected component is more dense 
than the bulk liquid ( 6). 

In any of the above mechanisms for rad ia ] compositional 
variations, the shape and composition of the solid-liquid inter- 
face is an important factor. Knowledge of the shape of the Inter- 
face at various points in the crystal during growth is required 
in order to determine the causes of radial compositional variations 
and to design growth procedures that minimize these variations. 

The high mercury vapor pressure in HgCdTe makes instrumentation of 
the ampoule very difficult; thus traditional means of interface 
demarcation by Peltier pulsing cannot be used. The purpose of this 
paper is to indicate a method by which the interface shape and com- 
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position may be inferred from estimates of the thermal field 
and measured composition data taken from the grown crystal. 

The consequences of interface curvature in a system with a 
highly sloping liquidus and solidus curve such as HgCdTe up- 
on compositional variations is indicated and experimental 
examples are presented. 

2 . 0 INTERFACE DECONVOLUTION - THEORY 

We wish to formulate the general problem of an alloy being 
directionally solidified at a constant pull rate R. Choose a 
co-ordinate system x,y,z fixed to the sample with the y-axis 
aligned with the sample axis and ant i-paralled with the gravity 
vector g (positive y is up, toward the liquid). The following 
equations relate the fluid velocity u and concentration C. 

V • u = 0 , ( 1 ) 

|f + ( u • 7 ) u = -l/p o VP+v V 2 u+g (p/p -1), (2) 

If + G-VC = DV 2 C, (3) 

where P is a pressure ,v the kinematic viscosityvp the fluid 
density, D the liquid diffusion constant and t is time. We 
have assumed that diffusion in the solid is negligible and 
that the material properties are independent of concentration 
and temperature except for the difference between liquid and 
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solid. The above equations represent the simplified Ober bee k- 
Bou.ssine.sq app roximat ion . 

The thermal field is best written initially in a furnace 
fixed frame x,y,z whose origin is at the furnace midplane. If 
V is the media velocity seen from the furnace frame, we have 

+ V-V T = k V 2 T , (4) 

with < the thermal diffusivity and T the temperature. Boundary 

( ; 

conditions for eqn . (4) must specify an initial value for the 

thermal field T(x,y,z,0), and either a temperature over some 
surface or the thermal flux across that surface. The thermal 
problem and its means of solution is well known, however pro- 
blems in solidification are complicated by the presence of the 
solid-liquid interface in the sample. At the interface we have 
an additional boundary condition 

k L ( V T • n ) L = k s (VT-n) s ' (5) 

with n the normal to the interface, k_ and k the thermal con- 

L s 

ductivity of the liquid and solid, respectively. Here we have 
assumed that the solidification velocity is sufficiently slow 
that latent heat effects are ignor able . The thermal field equa- 
tions may be transformed to the sample fixed co-ordinate system 
with x-x , y-y+Rt-y 0 , z-z , where v is the position of the sample 
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tip at t=0. We relate the media velocity V to the pull rate and 
fluid velocity u by 

V « R + u (6) 

where R ~(0,-R,0) and u=0 in the solid. The boundary conditions 
specifying the thermal field must be transformed into the sample- 
fixed coordinate system and in general will be a function of time. 
The resulting equation for the thermal field in the sample fixed 
co-ordinate system be c ome s 

~ + u- VT = k V 2 T , (7) 

and the boundary conditions must be transformed into the sample 
co-ordinate system. 

The boundary conditions for eqns. (1-3) may be written in the 
sample co-ordinate system as 

C(y^oo) = C q (8) 

v*n (C -C ) = p T /p D VC • n at the interface , (9) 

S i i u s 

C T = k C- at the interface , (10) 

s I I 
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at the interface. 


( 1 i ) 


v • n ( p -p ) = (u * n ) p 
L S L 

T ] . = f ( C £ ) (12) 

P= p(c,T) - (13) 

Equation (8) assumes an infinite sample with the bulk liquid 

solute concentration C . Equation (9) is the solute conservation 

o 

relation applied at the interface where v is the solidification 

velocity, C ^ denotes the solute concentration in the liquid at 

the interface, C T the solid solute concentration at the interface 

s I 

and the solid and liquid densities P and P T are also evaluated 

s L 

at the interface. Equation (10) relates C g j. to C^. via the equi- 
librium segregation coefficient k . Equation (11) is a condition 
on the normal component of the fluid velocity at the interface, 
and we also assume the tangential component of u vanishes at 
the interface (no -slip condition). The temperature at the inter- 
face is determined via the phase diagram, and equation (12) relates 
the liquid interfacial composition to the liquid us temperature 

via the function f, i.e. f is the 1 iq nidus curve of the phase dia- 
gram. Equation (13) expresses the density as some proscribed function 
of the concentration and temperature. We have assumed there is no 
undercooling and that capillarity effects are ignorabie. 

The set of equation (1-4) together with their boundary conditions 
represent the general description of an alloy undergoing directional 
solidification. Several authors (l,2) have addressed the solution 
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to these equations using per turbut ion techniques. In a later 
paper, we will approach the problem of shape instability utiliz- 
ing these equations. For the purposes of this paper, we wish to 
describe a different approach that allows, under certain circum- 
stances, the problem described above to be solved for the interface 
shape and concentration, without solving explicitly the rather 
formidable set of equations. 

We begin by assuming that we have analyzed a crystal grown 

by directional solidification and determined the solid solute 

concentration profile. For the sake of description,, we will 

assume that the crystal is azimuthally symmetric so that we have 

a two dimensional problem. Let C g -.(x,y) be the solid composition 

profile measured on the crystal. In actual practice, we measure 

the solid solute concentration at discrete points; a polynomial 

fit may then be used to construct c (x.y). 

s 

Futhermore, we assume that the furnace is sufficiently 
well characterized that a thermal model based on equation (4) will 
yield a thermal field sufficiently accurate for our purposes. At 
some point in time, we assume as an initial guess the position 
and shape of the solid-liquid interface. A thermal model together 
with appropriate boundary conditions is used to calculate a trial 
thermal field T(x,v), utilizing the assumed interface shape and 
position to apply the thermal boundary condition at the interface. 

We also use the solidus curve of the phase diagram S ( C g ( x , v ) ) , 
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where S is the solidus temperature for the solute composition 

C s (x , y ) and we assume no undercooling is present. We thus have 

two temperature distributions at any point in time in the crystal; 

the temperature at which each point froze S(C (x,y)) 5 and the 

s 

calculated temperature in the crystal at each point, T(x,y). 

The solid-liquid interface is the locus of all points whose 

solidification temperature S(C (x,y)) is equal to the temperature 

s 

T (x , y ) at that point in time. We have implicitly assumed that 
both S and T are monotonic functions, i.e. growth conditions are 

I 

such that there is only one interface and the liquid is never 
constitutionally supercooled. 

We may form a function 1= S ( C (x , y ) ) -T ( x , y ) which gives the 
k sV 

interface location x ,y by setting 

I(x*,y*) - 0, (14) 

that is, the zeros of I(x,y) are the interface co-ordinates at 
that point in time. 

The interface position calculated by this ’’deconvolution" 
procedure may not be the same as assumed initially to calculate 
the thermal field T(x,y). We now use the new interface co-ordinates 
in the thermal model to calculate a new thermal field. We should 
note that we now have an interface position and the temperatures 
along that interface; this is a boundary condition in calculating 
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the new thermal field. The new thermal field is used with the 

measured concentration data (which does not change) in the de- 

convolution procedure described above to calculate new interface 
* * 

co-ordinates x ,y , This iteration procedure continues until we 
find interface co-ordinates that do not change from the previous 
iterations. The final result is an interface location together 
with a self-consistent thermal field that are both consistent 
with the measured solid solute profile. The flow chart in Figure 
1 indicates the steps involved in the iteration procedure. 

f "■ 

In essence, this procedure substitutes experimentally mea- 
sured concentration data for solving the coupled solute equations. 

The thermal model must assume that heat transport due to fluid 
motion is ignorable. The iteration procedure retains the coupling 
of the solute equation and the thermal field which involves equations 
(3, 4, 9, and 12). Upon convergence, we have interface co-ordinates 
that are self-consistent with the boundary conditions for the thermal 
field and the experimentally measured solute profiles. We have not, 
however, determined the solute profile in the liquid ahead of the 
interface and thus have not solved for u in equations (1 and 2). 

Convergence limits of this iterative procedure have not yet 
been established. It is in the nature of a compact transformation, 
so that convergence will be assured within a region sufficiently 
close to the final valve. To date, this iterative procedure has 
not been carried out and thus the constraints on the initial guess 
for the interface location to insure convergence have not been de- 
termined. Examples of the deconvolution process and how various 
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thermal fields affect the interface shape in HgCdTe are given in 
the following section. 

3.° INTERFACE DECONVOLTION - EXAMPLE 

Although a complete treatment of the solidification problem 
would require carrying out the iteration process described above, 
we may illustrate the deconvolution method by a simple example. 

Figure 2 is a plot of measured isoconcentration lines converted 
into so lidif icat ion temperatures using the pseudo -binary phase 

■f " 

diagram for HgCdTe for a crystal grown under directional solidifi- 
cation. The solute concentration values were measured at discrete 
points along the crystal using an SEM with an energy dispersive 
spectrometer. This discrete data was then fit to a polynomial at 
each radial (x-value) point and interpolated to obtain isoconcen- 
tration lines. The concentration was converted via the phase diagram 
to an equivalent solidus temperature which is indicated on the dia- 
gram. These data are far from ideal, but do indicate the magnitude 
of radial segregation present in directionally solidified HgCdTe 
crystals. The abrupt drop in the iso co ncen t r a t ion lines at the right 
hand edge of the crystal cannot be explained in terms of the simple 
model presented here, we have ignored these data points in our 
analy s is . 

The central issue is that if we assume that the solid-liquid 
interface is of uniform composition and hence an isotherm we must 
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interpret the composition data in Figure 2 as indicating a 
extremely curved interface shape and consequently highly curved 
isotherms* However, there is no reason to assume that the inter- 
face is of uniform composition. Convection, lateral diffusion, 
and shape instability all may act to create an interface that is 
neither an isotherm nor an isoconcentrate* We may illustrate the 
magnitude of this effect by computing an interface shape consistent 
with the data in Figure 2 for various trial thermal fields. 

Figure 3 is a schematic of two typical thermal profiles that 
might be obtained in a crystal growth experiment. These thermal 
fields are not specific to any furnace or sample configuration but 
are used as examples only. They are, however, representative of 
the types of isother d is s t r ib u t ions obtained in HgCdTe from thermal 
models . 

The interface location is deconvoluted from the data contained 
in Figures 2 and 3. At selected radial (x-value) positions the 
axial composition data is converted via the phase diagram to a 
solidus temperature and fit to a least squares polynomial giving 
a function ( C (y ) ) for each selected x-value. In a similar manner 
the thermal data is fit to a polynomial yielding T (y) . We note 
that in general the thermal field expressed in the sample is a func- 
tion of time and sample pull rate. In this example we assume that 
the thermal field has attained its steady state value and that growth 
velocities are sufficiently slow that heat transfer due to sample 
motion does not perturb the thermal field. 
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Given S^(C ( y ) ) and T (y) we find the zero of X „ ( y ) - 

X S X X 

S x (C s (v)) -T_(y) for each x value which gives the interface 
* * 

co-ordinates x , y . The results are illustrated in Figure 4. 

For each assumed thermal field two different interface positions 
are shown; the number at the right hand side is the approximate 
location of the Interface in mm along the solidified sample. The 
mid-plane is the location of the flat isotherm in each thermal 
profile. The two interfaces do not exist at the same point in 
time of course, but rather are calculated at different times in the 
growth. The time interval between the interface at 32mm mud the 
interface at 101mm is 90 hours. The sample translation rate was 
0.762 mm /hr . 

Figure 4 indicates that despite the highly curved isocon- 
centration lines shown in Figure 2, we find that only a slight 
curvature in the interface results from the typical thermal fields 
assumed here. The implications are that very small interface curv- 
atures lead to gross radial segregation in HgCdTe. The first thermal 
field yields interfaces in a region of convex isotherms while the 
second thermal field gives an interface located in a region of 
concave isotherms. We note however that the interface is concave 
in each case due to the highly concave nature of the isoconcentration 
lines. In each thermal field assumed, we find that the interface 
moves downward from the initial transient region at 32mm along the 
crystal to a near steady state position at a distance solidified 
of 10 1mm. This is consistent with results obtained in a one d linen - 
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sional diffusion model ( 6). 


4.0 DISCUSSION 

We have described a technique by which the solid-liquid 
interface location and composition may be determined from 
compositional analysis of a crystal and determination of the 
thermal profile via suitable thermal models. The resulting 
Interface is consistent with both the experimentally determined 

f 

solute composition profile and the calculated thermal field. 

The technique allows for an interface that is neither an iso- 
therm nor an isoconcentrate, and can be used to determine 
the interface location at any point in the crystal. 

In order to be useful in determining the causes of radial 
segregation in systems such as HgCdTe , we must first ensure 
that the thermal model adequately describes the thermal con- 
ditions encountered during an actual experiment. Thermophysical 
properties of solid solution semiconductors such as HgCdTe are 
not well known, and the high vapor pressures encountered prevent 
a direct measure of the thermal profile with thermocouples. The 
iteration procedure described in this paper depends on accurate 
values for both the compositional profile and the thermal profile. 
The problems encountered in determining the thermal field are the 
same, however, as encountered in attempting to solve the set of 
coupled equations described in section 2. 
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Due to a lack of thermophysical properties for HgCdTe we 


have not carried out the full iterative procedure necessary to 
determine an interface consistent with the compositional profile 
and experimentally measured thermal boundary conditions. By way 
of example we have illustrated that thermal fields typical of 
those encountered in real growth situations leads to interface 
shapes with only a slight degree of curvature being consistent 
with measured iso concentration profiles that are highly curved. 
This indicates that in a system such as HgCdTe, mechanisms that 

i 

lead to even a slightly curved solid-liquid interface may result 
in gross radial compositional variations. One .-ueh mechanism, 
that of shape instability, will be addressed in a later paper. 
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FIGURE CAPTIONS 


FIGURE 1 

FIGURE 2 

FIGURE 3 
FIGURE 4 


Schematic Diagram of Interface Determination 
P rocedure . 

Solidus determined from Crystal Composition 
Data in HgCdTe . 

Thermal Profiles Applicable to HgCdTe Growth. 

Interface Locations in HgCdTe for Two Assumed 
Thermal Fields . 
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FIGURE 1. SCHEMATIC DIAGRAM OF INTERFACE DETERMINATION PROCEDURE 




